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We show that the distinct differences of low-lying THz-frequency dynamics between 
type-I clathrates with on-center and off-center guest ions naturally follow from a the- 
oretical model taking into account essential features of the dynamics of rattling guest 
ions. Our model analysis demonstrates the drastic change from the conventional dy- 
namics shown by on-center systems to the peculiar dynamics of off-center systems in 
a unified manner. We claim that glass-like plateau thermal conductivities observed for 
off-center systems stem from the flattening of acoustic phonon dispersion in the regime 
|fc| < |G|/4. The mechanism is applicable to other systems such as glasses or relaxers. 
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1. Introduction 

Type-I clathrate compounds with the chemical formula II 8 III 16 IV3o have recently at- 
tracted a great deal of attention in connection with the search of efficient thermoelectric 
materials. 1 ' This is because, though these compounds are crystalline, they exhibit glass- 
like phonon-thermal conductivities. 2 ~ 6 * ) Diffraction experiments 7-9 ^ and extended x-ray 
absorption fine-structure (EXAFS) studies 10,11 * 1 on these clathrates have revealed that 
rattling guest ions in cages take either the on-center or the off-center positions depend- 
ing on the relative size of the cage to the guest ion. It is remarkable that clathrates with 
"off-center" rattling guest ions exhibit glass-like low-lying modes in the THz-frequency 
region, 12 ^ which have been observed in specific heat measurements, 6 ' 13 - 1 inelastic neutron 
scattering experiments, 14,15 - 1 and optic spectroscopy. 16 ~ 20 ) 

There are several theoretical works on the THz-frequency dynamics of type-I 
clathrates. 21 ~ 23 ) To complement these works, we consider in this paper another physi- 
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cal model to investigate the differences between the on-center system (ONS) and the 
off-center system (OFS). Our approach emphasizes the difference of ONS and OFS on 
the emergence of peculiar low-lying THz-frequency modes observed in type-I clathrates 
with off-center rattling guest ions. In addition, we discuss the effect of the long-range 
electrostatic interaction between the rattling guest ions, which has been overlooked in 
other work. Recent experimental data in terms of far infrared spectroscopy at low tem- 
peratures by Mori et a/. 20 -* for OFS provide important information on the difference 
about THz frequency dynamics from ONS. 

We construct the equations of motion to determine the dispersion relations of acous- 
tic and low-lying optic modes for ONS and OFS of clathrates. The calculated results 
quantitatively clarify the distinct differences of the THz-frequency dynamics between 
ONS and OFS in the following points: dispersion relations of acoustic phonons, be- 
haviors of low-lying optic modes and characteristics of thermal transports driven by 
acoustic phonons. 

2. Description of the Model 

We illustrate compounds with the chemical formula II 8 III 16 IV3o as a prototype of 
type-I clathrates. Rattling guest ions in these compounds take the on- or off-center 
position in cages of 14 hedron depending on the ratio of cage size and ionic radius of 
rattling guest ion. For example, the rattling guest ions in 14 hedron in (3-Ba. 8 GsiiQSn 30 (j3- 
BGS) takes the off-center position at Uq = 0.43 A from the center of the tetradecahedral 
cage. 6 ) The deviation Uq = 0.43 A is 7.4% of the nearest-neighbor distance d = 5.84 A 
between Ba 2+ ions. The shape of cages restricts the motion of the off-center guest ions 
into quasi two-dimensional plane parallel to the hexagonal faces of the cage. Hence, the 
potential energy for off-center guest ions is shaped like the bottom of a wine bottle with 
fourfold hollows. 4 ' 24 ) On the other hand, the rattling guest ions in n-type /3-Ba 8 Gai 6 Ge 3 o 
(/3-BGG) take almost on-center positions in cages 4,9 ) forming a potential shaped like 
the bottom of a "wine glass" . 

We focus our attention on the lowest two bands (acoustic and optic modes) in the 
THz-frequency range. This is because these two modes are crucial for elucidating the 
mechanism controlling the differences. We treat the cage as having mass M and an 
effective charge e* c ] and the rattling guest ion, mass m and charge e* G . We express the 
position of the £-th cage at time t as Ri + Re{t), where Re is the equilibrium position 
of the £-th cage center. The vector rg(t) represents a small displacement from rg. The 
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position of the rattling guest ion from re is defined by the vector U e + Ue(t), where 
Ui is the equilibrium position of the rattling guest ion from re, and Ui(t) is a small 
displacement from Ue at time t. Note that Ue ^ in the case of off-center guest ion, 
while the on-center case becomes U e = 0. 

Scaffold cages compose cubic lattice, 7 ~ 10 ) which is invariant under translation by 
any lattice vector. Each cage is elastically connected with nearest neighbor cages. We 
denote the coupling strength by the harmonic force constants f» and f±, which repre- 
sent longitudinal (dilation) and transverse (shear) modes. With these quantities, the 
potential energy is expressed by 

f 2 



t,i 



fi. 

re,i - fe+i,i\ 



= 2^/ 7 |Q fc7 | 2 sin 2 (^), (1) 

k,i ^ ' 

where a is a lattice vector, and 7 denotes the species of three modes: ||, _L, and _!_'. The 
second line is obtained by the Fourier transformation: 



L£^*"'- (2) 



The rattling guest ions are trapped in anharmonic potentials, for which we keep up 
to the 4th order of small displacement in the expansion. This choice enables us to treat 
both of ONS and OFS in a unified manner. The anharmonic potential is given, in terms 
of relative displacement we = Ue(t) — re(t), by 



^ 2 = E 



^\U £ + we(t)\ 2 + ^\Ue + we(t)\' 



(3) 



Note here that the anharmonic potential is expressed by the relative coordinate wg. 
This allows us to legitimately treat the THz frequency dynamics of a coupled system 
consisting of guest ions and networked cages beyond the adiabatic approximation. In 
addition, it is possible to make the Fourier transformation for We on the fc-space since 
We is independent of Ue- 

The potential V2 involves two types of anharmonic potentials depending on the 
values of parameters £ and rj. Depending on positive or negative £, Eq. (3) represents 
wine-glass or wine-bottle type potential where rj is positive constant. In addition to 
V2, we should take account of the electrostatic interaction energy V3 between off-center 
ions; we shall discuss this later. 
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3. Potential for Off-Centered Guest Ions 

So far, our discussion has been quite general. We now restrict attention to the 
case of OFS, for which we impose the condition £ < and 77 > in Eq. (3) for 
V 2 . This parameter set yields a wine-bottle type potential, which takes its minimum 
^2,mm = ~ £ 2 / (4?7) at I U ^ = —£/r). The non-zero U e allows us to rewrite the potential 

(3) as 

v 2 = ^j2\ u t +w t&\ 2 

Here we define the fluctuation around the equilibrium position of off-center guest ions 
by introducing complex number given as 

We = eWt+M&y^) 

where We, is a complex number representing the vector U e + We- From Eq. (5), the 
deviation from the center of cage becomes [/oe j6 V The equilibrium angles {9e} are ori- 
ented at random reflecting glass-like behaviors observed in specific heats and thermal 
conductivities for OFS. 2 ~ 6 ) 

The form of Eq. (5) is analogous to that introduced in the local gauge theory in 
relevance with the Higgs mechanism, in which a covariant derivative plays a role for 
providing with masses to dynamical variables he(t) and 59e(t). We take a different 
approach to give physical implications for these dynamical variable as explained below. 

The random orientation of the vector U e. originates from four hollows of the anhar- 
monic potential 4, 24 ** in addition to inter-site dipole-dipole interactions between off-center 
guest ions. 25 ' 26 ) It has been confirmed 4 -' that off-centered guest ions experience a hinder- 
ing potential Vh(9) with a fourfold inversion symmetry along the azimuthal direction. 
The first-principles calculations have shown 14 to be ~ 20 K for Sr 2+ guest ions in 
/3-SGG. 24 ) 

The effect of long-range electrostatic force has been treated in our papers 25 ' 2&> for 
the elucidation of glass-like behaviours of OFS at low temperatures. Electric dipole 
moments of guest ions provide long-range interactions with dipoles in other cages. To 
make our argument clearer, let us illustrate as an example the two coupled dipoles 
p 1 and p 2 , whose potential is given by V12 = Vh{9i) + Vh{9 2 ) + Vi 2 (9i,9 2 ) with the 
dipole interaction V12 = (gqIIq) 2 cos(9i — 9 2 )/Rl 2 with the distance R\ 2 between the 
dipoles at 1 and 2. Two global minima (maxima) in (9i,9 2 ) configuration space appear 
at \9 l -9 2 \ = n (2tt). 



Ue + w e (t)\' 
2\UA 2 



(4) 



U + 



he{t) 
V2 J 



(5) 
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Fig. 1. Illustration showing the emergence of two freedoms of motion; stretching hi(t) and libration 

u 6e e (t). 



This argument can be straightforwardly extended to the case of multiple pairs pro- 
viding many potential minima in the configuration space V = (9\, 82, 63, ■ ■ ■ ), where 
the potential function is V123... = Yl^h(9i) + Yl^ij- This configuration acts as a new 
hindering potential in addition to the four- fold inversion symmetric potential Vh- The 
key yielding glass-like behaviors is that the dipoles p £ = e* G Ug in OFS constitute an 
equilateral-triangle structure among next-nearest neighbours dipoles and easily gener- 
ate a frustrated situation necessary to glass-like properties. The energy-scale of dipole- 
dipole interaction is estimated to be of the order of a few 10 K in Ref. 25 * 1 which is of 
the same order of the barrier height of Vh- 

This implies, in the temperature region ksT < Vh ~ a few 10 K, that off-center 
guest-ions execute angular fluctuation {59e} from the bottom of the potential at {9e}- 
Thus, We treat the phase factor as a dynamical variable trapped in a hindering potential. 
The variable he(t) in Eq. (5) represents small fluctuation along the radial direction 
directed at angle Qg. Thus, we have a small displacement of rattling guest ion at £-th 
site to the first order 

Wi (t) = ^=[hi{t) + iU 59 t (t)}e ie *. (6) 
v2 

We have depicted the situation in Fig. 1. 

The equation (4) becomes, under the situation | Ue\ > \ wg\ using Eq. (5) 

^ = ^£(W + &W)' (7) 

e 

with the definition of the effective force constants for stretching and libration modes 
given below 

l s = -2Ul + (h 2 £ )/^), & = -2&(l + (<W?)/4). (8) 
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The effective force constants £ s and £g in Eq. (7) are introduced to redefine £ in Eq. (4) 
in order to express the differences between angular motion and stretching motion. It 
should be emphasized here that the effect of random orientation {Ue} on effective force 
constants shall be involved in newly introduced force constant £ 7 in Eq. (11). 

We have employed the thermal average in these expressions. Thermal-averaged 
squared-displacements, (hf) and (58]), are proportional to temperature T. The Ra- 
man scattering 16 ' 17,19 ) and far infrared spectroscopy 18 ' 20 ) experiments show that spectra 
monotonically increase with temperature T. These observations indicate the relevance 
of anharmonic potential via (hj) and (S6j) since £ s and are directly related with the 
spectra as will be given in Eq. (19). 

It should be emphasized that anharmonic terms in Eq. (7) involved in £ s and £g 
become irrelevant below a few 10 K. Thus, we do not take into account the effect 
of anharmonic terms in the low temperature region discussed in this paper. In this 
connection, we remark that, at sufficiently high temperatures above a few 10K, guest 
ions execute free-circular motion and should recover the translational symmetry. 

4. Equations of Motion 

4-1 Coupling between Acoustic Phonons and Guest Ions for Off- Center System 

Let us discuss the effects of the random orientation of the equilibrium angles {6^} 
of off-centered guest-ions at the site I. There are two aspects on the effects to the THz 
frequency dynamics. 

The first is the coupling mechanisms between acoustic phonon modes and off- 
centered guest-ions trapped in cages. There are two possible types of the coupling: 
the symmetry-sensitive coupling, which is related to deformation of cages, and the 
symmetry-insensitive type related to the displacement of cages and guest ions. 

The potential function Eq. (4) for off-centered guest ions manifests the local po- 
tential interacting with cages. The situation is analogous to localized electron-phonon 
coupling in insulators. From the arguments of symmetry of modes for this case, Raman- 
active stretching mode couples with longitudinal acoustic phonons accompanied with 
"dilatation" , and the infrared active libration mode does with transverse acoustic modes 
without dilatation. Thus, when the symmetry-sensitive coupling is valid, the infrared 
active libration motion Uo89g primarily concerns with the transverse modes (_L) without 
dilatation, while the radial stretching hi(t) does with the longitudinal mode (||) with 
dilatation. 
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Fig. 2. Illustration showing the symmetry-insensitive coupling between acoustic phonons and guest 
ions. The figure (a) is the case that the off-center guest ion takes the far-side position for incident 
longitudinal acoustic phonons with the wave vector and (b) shows the guest ion sits perpendic- 
ular position to the same longitudinal acoustic phonons. The coupling becomes effective to parallel 
components of relevant displacements regardless of the transverse or longitudinal modes. 



The symmetry-insensitive coupling is independent of the symmetry of modes. This 
coupling becomes relevant to parallel components of each displacement regardless of 
the transverse or longitudinal modes. The situations illustrated in Fig. 2(a) and (b) 
provide the implication of this coupling mechanism. For example, consider the case 
that longitudinal acoustic phonons with polarization vector parallel to the y-axis are 
incident along y-axis when the guest ion takes the far-side position as depicted in 
Fig. 2(a). These longitudinal acoustic phonons dominantly couple with the stretching 
motion he(t) parallel to the polarization vector. Provided that the guest ion sits the 
right side along x-axis as shown in Fig. 2(b), the longitudinal acoustic phonons incident 
along y-axis couple with the libration motion Uq59i which is parallel to the polarization 
vector of the acoustic phonons. Thus, the coupling constant continuously spans from the 
libration one £g to the stretching one £ s . The same arguments hold for the coupling of 
transverse acoustic phonons. These cause the broadening of the spectral width through 
the distribution of force constants £9 < £ < £ s . 

Inelastic neutron scattering experiments are crucial to determine which coupling 
is dominant. We describe this point in detail when showing the calculated phonon- 
dispersion relations in Sec. V. 

The second is on the line-width broadening of the spectra of optic modes. Re- 
cent experimental data in terms of far infrared spectroscopy at low temperatures 20 ^ for 
OFS provide interesting features on the line-width broadening. They discovered that 
the spectral width of about 0.57 THz of the lowest-lying infrared active optic modes 
0.71 THz at 7 K decreases with increasing temperature. This feature contradicts with 
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the assumption that the anharmonicity of the potential Eq. (3) is a key element for 
interpreting the origin of the line-width broadening. 

These experimental results 20 ^ suggest that the random configuration of U i yields the 
broadening of optical spectra at low temperatures. We claim that the estimated spectral 
width is given by A ~ znl/4ire r f 3 , where f is the mean distance from off-centered guest 
ion to the nearest one, and z is a number of the order of unity representing the effective 
coordination number dominated by a few dipoles which happen to be closest. 

From the discussion made above, we reform the potential function Eq. (4) as, under 
the situation \ wi\ « \Ui\, 

1 — 4|— — | + anharmonic terms 

= - ^^£ 7 |u;^ 7 | 2 + const. (9) 

It should be emphasized here that the newly defined effective force constant £ 7 involves 
the effect of random orientation of {Ug} in a manner as explained after Eq. (11). We 
make the Fourier transformation for Eq. (9) using the relation, 

^w^E^' ( 10 ) 

which leads, in the form applicable to both cases of the symmetry-sensitive and the 
symmetry-insensitive couplings, 

Here, for the symmetry-sensitive coupling, the force constant £ 7 should be understood 
as £|| ~* £s an d £± — >* £0 by corresponding stretching motion to longitudinal mode and 
librational motion to transverse modes given in Eq. (8). For the symmetry- insensitive 
coupling, the force constant £ 7 in Eq. (11) should be independent of 7 and distributed 
from £q to £ s (£0 < £ 7 < £ s ) as mentioned in the above paragraph. Thus, we have 
incorporated the effect of randomly oriented Ui into the strength distribution of the 
coupling constants £ 7 . 

The relative displacement from the center of the cage induce dynamic electric- 
dipoles for both of ONS and OFS as discussed by Fano. 27 ^ The inter-site dynamic 
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dipole-dipole interaction V3 is expressed as 

e *2 

Vs = 22 9P3~[^ ' Wl ' ~ ' R M>){WZ' ■ Ruf)\, (I 2 ) 

where the dipole p e is defined by p e = e G W£, and Rw is a unit vector in the direction 
to the vector ru> = f£ — ?>• The Fourier transformed expression for Eq.(12) under the 
random phase approximation is obtained as 



O 2 |2 c ^.r nf 1 - 3 (gfc7- R 



of J 



4 Z^i**rl" ^3 

+ c.c, (13) 

where g fc7 means the polarization vector of the mode 7, and c.c. indicates complex 
conjugate of the first term. 

With using the Lorentz sums, Eq. (13) can be expressed as 

v* = -f-Y,\**\ L v ( 14 ) 

fe 7 

where L 7 = (qr fc7 ■ fc) 2 — | for small fc of cubic lattice and a; 2 = 47ms(e*) 2 /m is the 
squared plasma frequency of rattling guest ions with mass m and charge e G in cgs 
esu units, ub = 1/a 3 is its number density charge. We should note that L\\ = 2/3 for 
longitudinal mode and L± = —1/3 for transverse modes, respectively. 28 ^ 

The total kinetic energy of the cages (M J2 e \rg , 7 | 2 /2) and rattling guest ions 

( m E^, 7 (^,7) 2 / 2 ) is S iven by 

K = — 2 — + 7r2J^ 7 l 

fe,7 fc,7 
772 v — ■> / • • * • * \ 

+ yLl^-^7+ ■ Q kl ) , (15) 

fe,7 

where the Fourier transformation — > Q k is made by using Eq. (10). 

The equations of motion for two variables q k ^, Q kl with 7 = |, _L are obtained from 
the Euler-Lagrange equation as 

™4 7 (<?fc 7 + Qfe 7 ) = (£7 + m ^l L i)(l^ n (16) 

(m + M)w 2 7 Q fc7 = -mw 2 7 g fc7 + 4/ 7 sin 2 <3fc 7 ( 17 ) 

Here, the above equations of motion is given of the form applicable to both cases 
of symmetry-sensitive and symmetry-insensitive couplings. It is straightforward, for 
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example, to extend these equations to the symmetry-insensitive case by replacing £ 7 — > £ 

(£i<z<is). 

4-2 On-Center System 

For ONS, we impose the condition £ > 0, 77 > in Eq. (3). This yields Ue(t) = 0. 
For this situation, we have the Fourier transformed V2 as 

V2 = l y £,£\q* r (t)\', (18) 

fe,7 

where £ = £ + n(\wi(t)\ 2 ) /2. The characteristics of the force constants £ 7 in Eq. (11) for 
OFS and £ in Eq. (18) for ONS are quite different. The crucial difference is that £ 7 for 
OFS involves both the librational freedom and the radial one (£ s ), while £ for ONS 
is isotropic between longitudinal and transverse modes. The electrostatic interaction 
V3 for ONS becomes the same as Eq. (14). The equations of motion for two variables 
g fc7 , Q fc7 are obtained by replacing the parameter as £ 7 — > £ in Eq. (16). 

5. Phonon Dispersion Relations for OFS and ONS 

We determine the force constants / 7 and £ 7 in Eqs. (16) and (17) from the ve- 
locities of sounds and the data of optic spectroscopy. The force constants h , f± can 
be obtained from the sound velocities v» = 3,369 m/sec and v± = 1,936 m/sec for 
/3-BGS (/3-Ba8Gai6Sn 30 ) of OFS 29 ^ using the relation t> 7 = a^J / 7 / (m + M). Thus, we 
have f\\ = 18.1 N/m and f± = 5.77 N/m, employing masses m = 136 u for Ba atom 
and M = 8.55m fro cage, and a = 11.68 A. The force constants f\\,f± for /3-BGG 
(/3-BasGai6Ge3o) of ONS can be obtained using the sound velocities v\\ = 4,096 m/sec 
and v± = 2, 795 m/sec as well. 15 ) These yield f\\ = 26 A N/m and f± = 12.3 N/m with 
using M = 7.04m and a = 10.78 A. 

The force constants £ 7 in Eq. (11) and £ in Eq. (18) can be determined from the 
data of optic spectroscopy since the spectra for optic modes provide the information 
at fc=0, namely, T-point. Taking k — > in Eq. (17) and combining it with Eq.(16), we 
have the squared eigenfrequency for optic modes given by 



£ 7 + muipL 



in 



< - I 1 + iJ • < 19 » 

where 7 denotes the longitudinal or transverse polarization of optic modes. Equation 
(19) involves the squared plasma frequency u 2 arising from the fluctuation of charged 
rattling guest ions. We can estimate the magnitude of the plasma frequency to be 
u p /2it ~ 0.09 THz using mass (m = 137 u) and the charge of Ba 2+ (e* G = 2e) by taking 
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the number density ns = 0.628 x 10 27 m 3 and the relative electric susceptibility 
e r /e « 10, 25) the same one as that of Si crystal. 

We can estimate, using the expression Eq. (19), the force constant £ 7 for /3-BGS of 
OFS employing the data of Raman scattering 16 ) and far infrared spectroscopy. 18 * 1 Mori 
et al. 18 - ) have recently observed infrared active spectra at 7 K with the lowest-lying peak 
at 0.71 THz with line- width broadening of 0.57 THz for /3-BGS of OFS by means of THz 
time-domain spectroscopy. The peculiar line broadening with decreasing temperature 
found by Mori et al. 18 - ) is interpreted as the effect of random orientation of dipoles as 
mentioned in Sec. IV. 

The spectrum of A lff stretching mode coupled with longitudinal mode is not observed 
for /3-BGS of OFS due to technical reasons . 16 - 1 For other type of OFS, e.g., /3-SGG, the 
spectra at 2 K are obtained at 48 cm -1 , and the spectrum at 36 cm -1 for /3-EGG of 
OFS for Ai g mode . 16 - 1 We employ these to estimate the eigenfrequency of A lfl mode of 
/3-BGS as Uq/2it = 30 cm -1 (= 0.9 THz), which we use to evaluate the force constant 
£|| relevant to longitudinal optic mode. 

Under these assignments, we estimate ^ = 7.32(1 ± 0.25) N/m and £j_ = 2.21(1 ± 
0.25) N/m, using the estimated plasma frequency u p /2tt = 0.09 THz for Eq. (19). The 
contribution to optic eigenfrequencies from this plasma frequency is only 10% at T- 
point. The force constant £ for n-type /3-BGG of ONS is obtained at 2 K in a same 
manner using the data of Raman scattering uq/2tt = 32 cm -1 (= 0.96 THz) 16 - ) as 
£ = 8.23 N/m. 

We show, at first, the calculated dispersion relations in Fig. 3 for /3-BGG of ONS 
using the force constants obtained above for Eqs. (16) and (17). The acoustic phonon 
dispersions in Fig. 3 for /3-BGG are flattened below wo 7 for both transverse and longi- 
tudinal modes. We emphasize that the eigenfrequencies of optic modes at T-point are 
degenerate. The experimental data on the dispersion relations in terms of coherent in- 
elastic neutron scattering are available for n-type /3-BGG of ONS. 14 ' 15 ) Our calculated 
results shown in Fig. 3 are good agreement with the inelastic neutron scattering data 
for /3-BGG of ONS. 14 ' 15 ) 

Next, we present the calculated dispersion relations for /3-BGS of OFS for both of 
symmetry-sensitive and symmetry- insensitive couplings in Fig. 4 (a) and (b). As de- 
scribed in Sec. IV, these two couplings provide qualitatively different line-width broad- 
enings. 

It is remarkable in Fig. 4 (b) that the spectral width for symmetry-insensitive cou- 
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Fig. 3. (Color online) Phonon dispersion curves along the [001] direction for /3-BGG of ONS. Dis- 
persion relations for longitudinal modes are plotted with solid lines and for transverse modes with 
dashed lines. Dotted linear lines from the origin represent the long-wave limit of acoustic dispersion 
relations. The degeneracy of optic modes at T-point is observed at 4.2 meV arising from the on-cetered 
symmetric potential. 



(a) (b) 




0.25 0.5 0.25 0.5 
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Fig. 4. (Color online) Phonon dispersion curves along the [001] direction for /3-BGS of OFS. The 
figure (a) represents the symmetry-sensitive coupling; (b), the symmetry-insensitive coupling. Band 
dispersion relations for longitudinal (transverse) modes are illustrated by the region bounded by solid 
(dashed) lines in both of (a) and (b) . Dotted linear lines from the origin represent the long- wave limit 
of acoustic dispersion relations. Note that optic modes observed at T-point are separated at 2.2 meV 
and 4.1 meV in the case (a), but not for the case (b). 



pling becomes broader and gapless in the region of avoided crossing attributing to the 
coupling between acoustic modes and guest ions. While, the spectral broadening of the 
symmetry-sensitive coupling shown in Fig. 4 (a) reflects the width due to dipole-dipole 
interaction between randomly oriented dipoles. In order to determine which coupling 
is dominant, the experimental information of inelastic neutron scattering (INS) is re- 
quired. 

Quite recently, INS measurements for both /3-BGS of OFS and /3-BGG of ONS were 
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performed in the temperature range from 5 K to 290 K using the cold neutron disk- 
chopper spectrometer AMATERAS installed in Materials and Life Science Experimental 
Facility (MLF), Japan Proton Accelerator Research Complex (J-PARC). 30 ^ Preliminary 
results for /3-BGS at 5 K indicate that the line-width broadening for the lowest band is 
about 0.5 THz around the main peak of 0.53 THz. There is no clear gap in the density 
of states between optic and acoustic branches, which is distinctly different from the 
results for /3-BGG. These suggest that both coupling mechanisms are working. 

6. Thermal conductivities for OFS 

The frequency width at the avoided crossing is obtained using Eqs. (16) and (17) as 

Su n = 

where the last relation holds for a small contribution from the plasma frequency. This 
relation is valid for in our argument since the contribution from the plasma frequency is 
only 10% of ujq± for /3-BGS. Equation (20) indicates that the frequency 5u c ^ at avoided 
crossing is governed by the quantities uq^ and the square root of mass ratio \Jm/M. 
The frequency 5u n for OFS is much smaller than the case for ONS due to the inequality 

^ < i 

The formula of phonon thermal conductivity K p h{T) is given by 

K ph (T) = - / C phiJ (uj)vy£ 7 duj : (21) 
6 7 Jo 

where C p h tTl v 7 , and £ 7 are the specific heat, group velocity, and mean-free path of 
acoustic phonons of the mode 7, respectively. The cut-off angular frequency u n of 
the integral in Eq. (21) represents the cross-over frequency from extended phonons to 
quasi-localized phonons of the mode 7. For both of ONS and OFS, u n corresponds to 
the frequency at the avoided crossing since the eigenmodes of acoustic branches at the 
avoided crossing are strongly hybridized with the optic modes originating from rattling 
guest ions. Thus, we need to discuss the contribution from extended acoustic phonons 
below u n . 

Figure 4 (a) for /3-BGG of ONS shows that the avoided crossings occur at \k\ > 
\G\/4, which allows the Umklapp process ki + ft2 = £3+ G for acoustic modes consisting 
of elastic vibrations of cages. While, the avoided crossing of transverse acoustic modes, 
mainly contributing to phonon heat transport, for /3-BGS of OFS occurs at \k\ ~ | G\/A 
as given in Fig. 4 (b), indicating that acoustic phonons carrying heat are limited to those 
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with the wave number \k\ < \G\/4. 

The cross-over frequency u n is obtained by subtracting the the frequency 5uj c1 given 
by Eq. (20) as u cl = u 0j — 5u cl . This leads to 



Since the peak frequency of thermal distribution at the temperature T for extended 
acoustic phonons becomes hjj p ~ "i.&ksT taking into account the Stefan shift, our 
analysis indicates that the relation holds between the lowest optic eigenfrequency Woj_ 
and the onset temperature T p of the plateau expressed as, from Eq. (22) 



This is a general rule applicable to other types of type-I clathrates with off-center 
rattling guest ions. 

7. Conclusions 

We have investigated the THz-frequency dynamics of type-I clathrates by highlight- 
ing the differences between the on- and off-center systems. Concerning the low-lying 
optic mode, we have shown the energy range of the transverse mode in /3-BGS of OFS 
is much lower than that in /3-BGG of ONS, whose origin is due to the librational motion 
of off-center rattling guest ions. 

As for the dispersion relations, our calculated results shown in Fig. 3 are good 
agreement with the experimental data for /3-BGG 14,15 -' in terms of inelastic neutron 
scattering. The dispersion relations for /3-BGS given in Fig. 4 (a) and (b) show a wide 
fiat region for acoustic phonon branch of transverse modes arising from the low-lying 
librational optic mode. About the heat transport, we have found that the distinct 
difference of thermal conductivities between /3-BGG and /3-BGS is attributed to the 
position of avoided crossing in the Brillouin zone, where the hybridization between 
optic modes and acoustic modes occurs. 

To conclude, our model analysis serves as gaining deeper insight into the physical 
origins of peculiar low-lying THz-frequency dynamics observed for type-I clathrates 
with rattling guest ions. The formula Eq. (23) between u^± and T p can used as a 
rule for designing efficient thermoelectric materials with suppressed glass-like thermal 
conductivities. In addition, we claim that the physical mechanism proposed in this 
paper is applicable for elucidating the origin of the boson peak and plateau thermal 
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conductivity observed in glasses 12 '' or relaxers. 31 ^ 
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